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Complete determination of the polarisation state of light requires at least four distinct projective 
measurements of the associated Stokes vector. Stability of state reconstruction, however, hinges 
on the condition number n of the corresponding instrument matrix. Optimisation of redundant 
measurement frames with an arbitrary number of analysis states, m, is considered in this Letter 
in the sense of minimisation of n. The minimum achievable n is analytically found and shown to 
be independent of m, except for m = 5 where this minimum is unachievable. Distribution of the 
optimal analysis states over the Poincare sphere is found to be described by spherical 2-designs, 
including the Platonic solids as special cases. Higher order polarisation properties also play a key 
role in nonlinear, stochastic and quantum processes. Optimal measurement schemes for nonlinear 
measurands of degree t are hence also considered and found to correspond to spherical 2t-designs, 
thereby constituting a generalisation of the concept of mutually unbiased bases. 
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Measurement of the polarisation of light is a common 
problem in many fields of physics including quantum 
information, astronomy, quantitative biology and single 
molecule orientational imaging [1-4]. Typically, determi¬ 
nation of the polarisation state of light, as parameterised 
by a 4 X 1 Stokes vector, S, follows by making projec¬ 
tive measurements onto a set of known analysis states, 
with complete state reconstruction requiring a minimum 
of four distinct measurements [5] . Arguably the simplest 
so-called complete polarimeter is that which comprises 
of linear polarisers, oriented at 0°, 45° and 90° to some 
reference axis, and a circular polariser. By virtue of the 
linear nature of the measurement process, the intensi¬ 
ties transmitted through each polarisation state analyser 
(PSA), denoted by D, can be related to the incident 
Stokes vector via D = AS, where A is known as the 
instrument matrix (in this example D is 4 x 1 and A 
is 4 X 4 in dimension). Although intuitively simple to 
understand and easy to implement, this polarimeter per¬ 
forms sub-optimally. Much effort has been invested over 
the years to optimise the geometry of polarimeters using 
metrics such as the total variance on the inferred Stokes 
vector [6, 7], information content [8-11], the determinant 
of the instrument matrix [12-14] and signal to noise ra¬ 
tio [15]. Perhaps most popular, however, is the condition 
number, n, of the instrument matrix [13-20] which de¬ 
scribes the stability of the polarisation inference problem 
regardless of reconstruction algorithm and bounds the ex¬ 
tent to which relative measurement errors are amplified 
during state reconstruction. Smaller condition numbers 
imply more robust measurements. Use of a measurement 
set of greater than four analysis states, however, is known 
to mitigate the effects of noise [8, 20, 21]. Nevertheless, 
only limited consideration has been given to optimisa¬ 
tion of these systems. Drawing from results in discrete 
computational geometry, linear algebra and state tomog¬ 
raphy this work hence considers the optimisation of such 


measurement schemes in the sense of minimisation of n 
and discusses the associated geometric interpretation. It 
is established analytically that the minimum condition 
number is independent of the number of analysis 
states. Formal equivalence between minimisation of k., 
maximisation of the determinant of the associated Gram 
matrix and minimisation of the equally weighted variance 
is established. Optimality constraints are further derived 
and used to construct some illustrative optimal measure¬ 
ment sets. More specifically, it is found that the distribu¬ 
tion of optimal analysis states over the Poincare sphere is 
intimately related to spherical 2-designs. Accordingly, in 
contrast to previous reports, optimal measurements are 
shown not to necessarily correspond to inscribed polyhe- 
dra of maximal volume. 

Going beyond linear reconstruction of the Stokes vec¬ 
tor, measurements of more complex functions of the 
Stokes vector, D(S), are made in a number of applica¬ 
tions. Nonlinear light scattering and material characteri¬ 
sation, for example, gives rise to intensities which depend 
on products of Stokes parameters [22, 23]. Furthermore, 
random media and rough surfaces can be studied through 
the changes in the statistical properties of the polari¬ 
sation of light induced upon transmission or reflection 
[24-26]. Full characterisation, however, requires deter¬ 
mination of the underlying probability distribution func¬ 
tion, or equivalently all higher order statistical moments 
of the Stokes parameters. Optimal analysis states for 
such higher order problems are also considered in this 
Letter and their relationship with spherical t-designs es¬ 
tablished. Importantly, although the language of clas¬ 
sical polarimetry will be used throughout, it should be 
noted that the results given are equally applicable in the 
quantum regime for states of a given number of photons. 
For example, knowledge of higher order moments of the 
Stokes operators can give insight into hidden “quantum” 
polarisation in a classically unpolarised state [27, 28]. 
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FIG. 1: (a)-(e) Optimal measurement frames for m = 4,6,8,12 and 20 defining the Platonic polyhedra inscribed in the 
Poincare sphere. Corresponding analysis states lie on the family of cones shown, (f)-(k) Example optimal measurement frames 
and associated non-Platonic polyhedra for m = 6, 8, 9, 7 and 11. 


The action of a PSA on incident light of arbitrary 
polarisation, can be considered as a projective measure¬ 
ment, whereby the output light has polarisation match¬ 
ing the nominal analysis state of the PSA (denoted 
A = (1, Ai, A 2 , ^ 3)^/2 with = 1) inten¬ 

sity of • S, where A has been normalised to ensure 
the PSA is passive [9]. Each row of the instrument ma¬ 
trix associated with m > 4 different measurements is 
given by Aj {j = 1, 2,..., m) and A thus has dimension 
m X 4. The corresponding vector of measured intensities, 
D, is hence m x 1 . Geometrically, each Aj can be con¬ 
sidered as defining a point on the surface of the Poincare 
sphere through the reduced vector aj = (Aji, Aj 2 , 
such that A defines the vertices of a polyhedron inscribed 
within the unit sphere. Moreover, the set of vectors {aj} 
constitutes a measurement frame [15]. 

The condition number of the instrument matrix is ex¬ 
plicitly defined as /^ = ||A||||A+||, where A+ denotes the 
generalised inverse of A and || • || denotes the matrix 
norm (taken as the Hilbert-Schmidt norm throughout 
this work). Noting that the normalisation imposed on 
the rows of A requires Aj • A j = 1 / 2 , it follows that 

||A|p = tr[A^A] = m/ 2 . ( 1 ) 

For any given experimental setup ||A|| is thus constant. 
Accordingly, minimisation of the condition number is 
equivalent to minimisation of ||A+|p, an alternative fig¬ 
ure of merit known as the equally weighted variance 
(EWV) [ 6 , 20]. The EWV quantifies the noise ampli¬ 
fication in the reconstructed Stokes vectors assuming a 
least norm reconstruction and equal magnitude errors 
on each measurement [ 6 ]. Similarly to above ||A+|p = 
tr[(A+)^A+] = tr[®“^], where ® = A^A is a 4 x 4 ma¬ 


trix and denotes its inverse. The inverse of ® exists 
if 3ij are not all co-planar and moreover can be written in 
the form = adj[®]/|®| where | • | and adj[-] denote the 
determinant and adjunct of a matrix respectively. The 
condition number (/^ > 0) can thus be expressed as 


= tr[®]tr[®“^] = 


1 m 



m din |®| 

2 dBii 

1=1 


( 2 ) 


where Jacobi’s formula has been used [29] and Bij is the 
(i,jf)th element of ®. Upon differentiation of Eq. (2), 
application of the product rule and back-substitution it 
can be shown that 2d\iiK = —dln|®|. Minimisation of 
the condition number of the instrument matrix is hence 
also equivalent to maximisation of the determinant of ®. 
Geometrically, it is interesting to note that |®| represents 
the volume squared of a 4-parallelotope in For the 
special case of m = 4, maximisation of |®| is equivalent 
to maximising the volume of the tetrahedron whose ver¬ 
tices in IZ^ are defined by (see Fig. 1(a)) as has been 
previously reported [16, 17]. 

Hadamard’s inequality [29] states that the determinant 
of ® is upper bounded by the product of its diagonal 
elements. The maximum determinant is thus obtained 
when ® is diagonal whereby the diagonal elements also 
correspond to the eigenvalues pi (/ = 1 ,... ,4). Explic¬ 
itly ® can be expressed in the form ® = A^Aj, 

whereby it follows by inspection that Bn = pi = m/4. 
Furthermore maximisation of |®| is subject to the con¬ 
straint tr[®] = m/2 (c.f. Eq. (1)). Use of the method of 
Lagrange multipliers then yields Pi/^ = P 2 = P 3 = Pa = 
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m/Vl, i.e. the condition number of A is minimised when 

/ 3 0 0 0 \ 

m 0 1 0 0 , . 

®-l2 0 0 1 0 • 

\ 0 0 0 1 / 

Substitution of Eq. (3) into Eq. ( 2 ) gives n = ~ 

4.47214. The condition number of an optimised po- 
larimeter is hence independent of the number of analysis 
states m. It also follows that |®| = m^/6912 and the 
EWV = 40/m. The fall off in the EWV with increas¬ 
ing number of measurements reflects the noise reduction 
arising from greater measurement redundancy. 

Generalisation of the above to the problem of re¬ 
construction of an A/" > 3 dimensional vector Sat = 
{So,Si,...,Sn-iV, for which Sf, using 

m > N projective measurements onto analysis states of 
the form An = ( 1 , ..., can also be easily 

performed. Consideration of the N = 3 (i.e. reduced 
dimensionality) case is applicable to linear polarimetry 
for example, whereas higher dimensional generalisations 
are relevant to polarimetry of three dimensional fields, 
whereby Stokes vectors become 9 x 1 in size [26, 30]. hol¬ 
lowing the steps given above in the N dimensional case 
it is found that |®| is maximised when ® is diagonal with 
non-zero elements of /3i = m/4 and Pi = Pi/{N — 1) for 
/ 7 ^ 1. In turn it follows that |®| = (m/4)^(A^ — 1)^“^, 
^2 ^2N^-4N^4 and the EWV = 2 /A:Vm. 

Whilst the above treatment has considered the mini¬ 
mum achievable condition number and EWV, the opti¬ 
mal measurement basis has not yet been determined. To 
this end, Eq. (3) must be invoked which upon generali¬ 
sation implies the set of polynomial constraints 

m m 

Y^aj = 0 and E (4) 

i=i i=i 

where I at is the N x N identity matrix. A measurement 
frame is optimal iff Eq. (4) is satisfied. When V = 3 it 
can be shown [31] that the optimal measurement frame 
corresponds to a.j defining a regular polygon inscribed in 
a unit circle. Incidentally, the regular inscribed polygons 
have the maximum area of all inscribed polygons. Eor 
the N = 4 case (which is exclusively considered hence¬ 
forth) , it is found that Eq. (4) is satisfied if the set of vec¬ 
tors {a^}, or equivalently the vertices of the underlying 
polyhedron, constitute a spherical 2 -design in and are 
thereby closely related to mutually unbiased bases [32]. 
Proof of this result follows from the definition of spheri¬ 
cal t-designs as a collection of m points on the surface of 
the unit sphere in IZ^ for which the (normalised) integral 
of any polynomial, ^(S), of degree t or less is equal to the 
average taken over the m points [33, 34]. Use of the poly¬ 
nomial functions g = Sj and SjSk {j and k = 1,2,3) in 
this definition yields Eq. (4) . It is also worthwhile noting 
that a spherical t-design is also a t — I design [33] . 


Numerical codes can be used for determination of 
spherical t-designs in general [35], however, in view of 
the symmetry of the N = 3 solution, an analogous sym¬ 
metry in the N = 4 case is expected and can be used to 
guide the construction of some simple 2 -designs. Specif¬ 
ically, letting m = rs, and adopting the initial ansatz 
a.j = (sin 6 >g cos[(/)p + 4>g], sin 6 >g sin[ 0 p + 4>g], cos 6 >g)^ for 
j = I,...,m, p = I,...,r, q = I,..., 5, where pp = 
271 pjr^ Eq. (4) reduces to 

s s 

'^^cos0q = 0 and 0q = -. (5) 

,=i ,=i ^ 

In general, the solution to Eq. (5), and hence the choice 
of optimal measurement frame, is not unique (even allow¬ 
ing for the intrinsic rotational freedom). Nevertheless, a 
number of solutions can be found as is now illustrated. 

Even m: Eor even m a simple optimal frame follows 
by taking s = 2 whereby Eq. (5) is trivially satisfied when 
cos^i = — cos ^2 = I/\/3- Without loss of generality 4>i 
can also be set to zero. If 4>2 = 0 it follows that, de¬ 
pending on the relative sign of s'mOi and sin 6 ^ 2 , optimal 
frames possessing either an inversion or mirror symmetry 
can be generated. In both cases, however, the vectors 
lie on two cones with apex at the origin and apex an¬ 
gle 20q (see e.g., Eigs. 1(b), (c), (f) and (g) for m = 6 
and 8 ). Geometrically, taking 4>2 7 ^ 0 corresponds to a 
rotation of the measurement states on the lower cone rel¬ 
ative to those on the upper cone as depicted in Eigs. 1(f) 
and (g). Solutions describing the regular tetrahedron 
(m, r, s) = (4, 2, 2), octahedron ( 6 ,3, 2) and cube ( 8 ,4, 2) 
can be generated in this manner (see Eigs. I(a)-(c)). 

Going further, a measurement frame containing m 
analysis states can be partitioned into M subsets of size 
rrii = 2p^, where m = 2 /i^ = 2 g. Optimal bases 

can then be constructed (4>g = 0 is assumed henceforth 
for simplicity) by constraining the vectors of each subset 
to lie on cones with coinciding axes, but with differing 
apex angles. Explicitly it can then be shown that the 
apex angles Oi {i = I,...,M) of each cone are given 
by sin^ = 2pXi/{3gi)^ where 1 [31]- Note 

that in the limiting case of ^ 0 (whereby necessarily 
Pi = I), the corresponding cone collapses to its axis (e.g. 
Eig. 1(d)). Both the regular icosahedron and dodecahe¬ 
dron can hence be constructed (Eigs. 1(d) and (e)). Pos¬ 
sible optimal measurement frames for the m = 4, 6 , 8 ,12 
and 20 cases are thus defined by the vertices of the Pla¬ 
tonic solids inscribed in the Poincare sphere, in agree¬ 
ment with the numerical results of [20]. The inscribed 
polyhedra generated for m = 4,6 and 12 correspond to 
the polyhedra of maximal volume. In contrast, however, 
noting that inscribed polyhedra of maximal volume are 
Euclidean simplexes [36], i.e. each face is triangular, the 
cube and dodecahedron do not have maximum volume 
in contrast to the claim of [20]. Non-uniqueness of the 
solution to Eq. (4) implies the associated polyhedra for 
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arbitrary m are also not of maximal volume in general. 

Odd m: When m is a factorable odd integer the ansatz 
used thus far is also capable of generating optimal mea¬ 
surement frames. In Fig. 1(h), for example, a possible op¬ 
timal frame is shown for (m, r, s) = (9, 3, 3), with cos O 2 = 
0 and cos ^ 1^3 = ±l/\/2. For prime m alternative solu¬ 
tions must, however, be sought. It is well known that no 
spherical 2 -design, and hence no optimal frame, exists for 
m = 5 [34] , as can be verified by calculation of the Groeb- 
ner basis [37] of the set of polynomial equations given 
by Eqs. (4) and the constraints aj • a^ = 1 . Solutions 
for larger prime m can nevertheless still be found. One 
such solution set (also valid for factorable odd m) takes 
the form a^ = (sin Oj cos 4>j , sin Oj sin (j)j , cos Oj ), where 
(j)j = 27rjl{m — 1 ), cos Oj = ci for j = 1 ,3,..., m — 2 , 
cosOj = C 2 for j = 2 ,4,..., m — 1 , cos^^ = 1 and 
Q = [3 ± {—iy^3m{m — 4)]/[3(l — m)] as follows from 
Eq. (4). This solution is depicted in Eigs. l(i)-(k) for 
m = 7,9 and 11. 

The underlying mathematical framework found above 
hints at a potential extension of the optimisation proce¬ 
dure to higher order measurement problems. Consider, 
for example, the measurement of a nonlinear function 
D{S), such as the generated intensity in an optical non¬ 
linear conversion problem. Letting P/c(S) denote a com¬ 
plete basis of polynomial functions, ordered by increasing 
polynomial degree and indexed by k = 0,1,... ,00, the 
measurand can be decomposed according to 

00 

D{S) = J2FkPk{S) (6) 

k=0 


where Fk are the associated expansion coefficients which 
are now the unknowns of interest. are assumed to be 
orthonormal over the Poincare sphere. Practically, a fi¬ 
nite number of measurements, m, are made in directions 
{A,} SO as to sample D(S). Eurthermore, the sum in 
Eq. ( 6 ) must be truncated at a finite order P, such that 
Eq. ( 6 ) reduces to D = PF, where F = (Fq, Pi, • • • 5 Pk)^ 
and [^]jk = Pk{-^j)‘ In general, at least K F 1 mea¬ 
surements are required for complete determination of F. 
Optimality in this case can again be quantified using the 
condition number of P. Assuming the polynomials 
satisfy P^(S)P/c(S) = C, where C is a constant, 

a trace constraint analogous to Eq. (1) follows. Noting 
Po(S) is also a constant, similar derivations to above yield 


_ C CK‘^ 
“ a + c - p2 


(7) 


where the corresponding optimal Gram matrix PiP must 
be diagonal. The resulting constraints (c.f. Eq. (4)) 
are satisfied if the analysis states constitute a spherical 
2 t-design, where t is the degree of the polynomial Pk- 
Proof follows in a similar fashion to above. As a con¬ 
crete example, if the polynomial functions Pk are taken 


as the spherical harmonics Yim (whereby the k index in 
Eq. ( 6 ) denotes a suitable lexicographic ordering of the 
indices up to maximum degree / = t, it follows 

that K = [t F 1)^ and Pq = l/V^. Using the addition 
theorem = (2/ + l)/( 47 r) the minimum 

achievable condition number is found to be /^ = (t + 1 )^. 
Importantly, it should be noted that spherical 2t designs 
do not necessarily exist for arbitrary m [35] . In general it 
is therefore found that estimation of higher order prop¬ 
erties not only becomes more ill-conditioned, but also 
requires a larger number of measurements to be made. 
Use of optimal measurement sets is therefore critical to 
reconstruction quality in this case. In this sense the Pla¬ 
tonic solids perform well as measurement frames since 
they constitute higher order t-designs, e.g. a regular do¬ 
decahedron is a spherical 5-design, whilst possessing a 
relatively small number of analysis states. 

In summary, optimal measurement frames for the re¬ 
construction of the Stokes vector S of polarised light have 
been analytically and geometrically investigated. This 
analysis can also be applied to the input polarisation 
states in Mueller matrix polarimetry [38] . Equivalence of 
optimisation based on the EWV, the condition number n 
of the associated instrument matrix and the determinant 
of its Gram matrix was established. Gonstraints on the 
optimal analysis states were derived and found to be sat¬ 
isfied by states defining spherical 2-designs. It followed 
that minimisation of k does not necessarily correspond 
to maximisation of the volume of the corresponding in¬ 
scribed polyhedron. Einally, results were extended to 
consider optimal frames when the measurand is a polyno¬ 
mial function of S of degree t. In this case optimal frames 
correspond to sets of analysis states constituting spheri¬ 
cal 2t designs. This work provides the means for optimal 
polarisation state tomography and hence paves the way 
for practical study of nonlinear or stochastic properties of 
polarisation, which are of interest in both biological and 
physical contexts. Within the former, for example, non¬ 
linear polarisation studies can provide insight into cellu¬ 
lar and molecular structure [23, 39], whilst the latter can 
enable study of fundamental quantum polarisation prop¬ 
erties [28, 40, 41]. Einally, these results also present in¬ 
teresting opportunities for establishing optimal schemes 
for quantum tomography of two level systems due to the 
underlying geometric parallels with the Bloch sphere. 

The authors acknowledge Markus Grassl for suggesting 
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for insightful discussions. AE acknowledges financial sup¬ 
port from the Gordon and Betty Moore Eoundation. 
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